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Abstract 

Self-magnetically  insulated,  high  voltage  transmis¬ 
sion  lines  are  used  in  inertial  confinement 
fusion  particle  accelerators^  to  transmit  power 
from  the  vacuum  insulator  to  the  diode.  Injection 
and  output  convoluted  sections  pose  special 
problems  in  establishing  the  desired  electron  flow 
pattern  needed  to  maintain  high  overall  efficiency. 
A  time  independent,  2-D  numerical  code  for  planar 
or  triplate  geometries  calculates  the  motion  of  a 
test  electron  through  the  tapered  input  or  output 
convolutes.  The  1-D  parapotential  model  is  assumed 
to  be  appropriate  at  each  position  and  the  magnetic 
field  and  potential  distribution  are  calculated 
in  the  vicinity  of  the  particle.  The  electric 
field  is  then  calculated  from  Gauss's  Law,  and 
the  electron  motion  is  calculated  relativistically . 
The  results  show  that  the  electron  canonical  momen¬ 
tum  in  the  direction  of  flow  changes  as  the  elec¬ 
tron  passes  through  a  convoluted  geometry.  As 
shown  by  Mendel,  these  electrons  flow  between  the 
conductors  after  the  convolute  without  re-inter- 
secting  the  cathode.  We  hypothesize  that  these 
electrons  lead  to  the  losses  observed  in  long 
self-magnetically  insulated  lines.  Results  of 
calculations  are  correlated  with  results  of  the 
Mite  power  flow  experiment. 

Introduction 

Transition  sections  into  the  magnetically  insulated 
transmission  lines  can  excite  an  apparent  instabi¬ 
lity  in  the  electron  flow  within  the  transport 
sectionJ,£t  and  cause  severe  energy  losses.  A 
numerical,  time  independent  2-D  code,  MITL,  has 
been  written  to  investigate  the  effects  of  these 
transition  sections  on  the  flow  pattern  within  the 
transport  section.  The  code  is  used  to  examine 
the  input  transition  in  the  Mite  experiment.  ’’ 

The  results  suggest  that  the  input  transitions  pro¬ 
duce  electron  flow  in  which  the  axial  canonical 
momentum  Px  is  approximately  lO"^  kg-m/s  or  10~° 
of  that  allowed  in  the  line.  The  transitions  that 
produce  broad  canonical  momentum  distributions 
F(P'X)  are  correlated  with  efficient  power  transport 
in  the  experiments.  Those  that  produce  narrow 
distributions  are  correlated  with  lossy  transport 
experimentally . 


*This  work  was  supported  by  the  U.S.  Dept,  of 
Energy,  under  Contract  DE-AC04-76-DP00789. 


Description  of  the  Problem  and  the  Approach 

The  rectangular  geometry  used  in  MITL  is  appro¬ 
priate  to  triplate  type  transition  sections  and  is 
i.  a  own  in  Fig.  1.  The  effective  width  w  of  the 
lines  perpendicular  to  E  and  E  x  B,  and  the  con¬ 
ductor  separation  d  are  variables  that  specify  the 
geometry  of  either  the  input  or  the  output  transi¬ 
tion  section.  Time  independence  is  assumed  in 
MITL  since  the  transit  time  through  the  transi¬ 
tion  sections  are  typically  a  few  nanoseconds  and 
are  much  shorter  than  the  pulse  duration.  Electro¬ 
magnetic  fields  within  the  magnetically  insulated 
lines  are  calculated  on  the  assumtion  that  the  1-D 
parapotential  equations  derived  by  Creedon0  are 
appropriate  for  each  position  in  the  2-D  convolute. 
The  justification  for  using  the  1-D  parapotential 
model  is  two  fold:  First,  the  scale  length  for 
the  geometry  variation  is  long  compared  to  the 
separation  between  conductors  so  the  flow  is 
approximately  one  dimensional  at  each  position. 
Secondly,  the  parapotential  theory  adequately 
describes  the  relationship  between  the  total  cur¬ 
rent  IT,  the  boundary  current  Ij  inside  the 
negative  conductor,  and  the  voltage  VQ  for  a  given 
line  with  a  vacuum  wave  Impedance  Z  .  The  agree¬ 
ment  with  2-D  electromagnetic  PIC  simulations0 
and  experiments0’6  is  within  the  numerical  and 
experimental  uncertainties  respectively.  The 
electromagnetic  fields  from  the  2-D  calculation0 
for  a  Mite  like  line  at  2.4  MV  with  Iy  =  450  kA 
and  Ig  =  243  kA  in  a  coaxial  geometry  with 
ZQ  =  8 fl ,  have  been  compared  in  Refs.  6  and  7,  and 
the  agreement  justifies  the  use  of  parapotential 
theory  in  these  calculations . 


POWER  FLOW  — — - ► 


Fig.  1.  Transition  section  geometry  in  MITL.  The 
line  width  is  specified  in  the  Z  direction 
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The  test  electron  is  then  injected  into  the  con¬ 
volute  and  its  motion  followed  through  the  convolute 
and  into  the  uniform  self-magnetically  insulated 
transmission  line.  Once  inside  the  uniform  line, 
the  canonical  momentum 

Px  =  7emUx  -  eA*  (1) 

is  a  constant  of  its  motion,  where  m  is  the  elec¬ 
tron  rest  mass,  -e  is  the  electron  charge,  Ux  is 
its  axial  velocity,  A^.  is  the  axial  component  of 
the  magnetic  vector  potential  and 

\  =1/(1  -  U2/c2)1/2  (2) 

O 

for  an  electron  of  speed  U  and  c  «  3  x  10  m/s. 

The  sum  of  its  kinetic  and  potential  energy  is 
also  a  constant  of  its  motion.  Since  the  problem 
is  assumed  to  be  electrostatic,  the  energy  is  not 
changed  by  the  convolute.  However,  since 
dL/dx  /  0  in  the  convolute,  then  Px  is  changed 
by  &PX  according  to  Lagrange's  equation  as  the 
electron  moves  through  the  convolute.  Conse¬ 
quently,  the  problem  is  reduced  to  calculating  aPx 
accurately.  Generally,  the  limitations  of  finite 
cell  size  and  a  finite  number  of  particles  in  2-D 
self-consistent  simulations  severely  limit  the 
accuracy  with  which  aPx  can  be  computed.  The 
numerical  noise  is  avoided,  at  the  expense  of_self- 
consistency,  by  using  analytic  equations  for  S  and 
8.  Without  self-consistency  the  calculations  are 
uot,  however,  quantitatively  exact.  The  value  of 
these  calculations  is  the  insight  they  provide 
into  the  effect  of  convolutes  on  the  electron 
flow. 


The  initial  value  of  X  =  XQ  is  chosen  for  each 
calculation. 


The  total  current  flowing  through  a  magnetically 
insulated  line  is  given  by 


V  +  (V  -  2>1/2]  +  ~T-_  y-~in] 

'  'm 


(4) 


y  =  —  2  +  i  (5) 

m0C 

where  is  the  value  of  Y  at  Ym  the  edge  of  the 
flow  pattern,  and  Y  is  given  by  the  applied 
line  voltage.  The  local  line  geometry  determines 
the  value  of 


60  w(x) 

8  =  =  ZzdUd  (l 

o 

The  input  parameters  are  the  line  profile  w(x)  and 
d(x)  for  each  axial  position  x,  the  voltage  VQ  at 
the  anode  for  V  =  0  at  the  cathode,  and  the  total 


current  I-  through  the  structure.  For  a  given 
position  tx,y) ,  the  value  of  y  (x)  is  calculated 
from  Eq .  4.  The  voltage  V(x,y)  is  given  by 

me2  Y 

V  -  (cosh  -g - 1)  0  <  Y  <  Yb  (7a) 


or 

me2  ,%<**■>  -  ^m  <Vd>  , 

v  -  —  (  — — r=T - )  Y~  <  Y  <  Y- 


m 


m 


where 

and 


ci  - 

Ym  -  Cj_  Zn  [Ym  +  (7m2  -  1)1/2] 


(7b) 

(8) 

(9) 


Description  of  the  Program 

The  program  involves  the  choice  of  the  initial  con¬ 
ditions  for  the  electron  as  it  leaves  the  cathode 
plasma  in_the  convolute, _the  calculation  of  the 
electric  E  and  magnetic  B  fields  in  the  vicinity 
of  the  test  electron,  and  the  integration  of  the 
relativistic  equation  of  motion  as  the  particle 
progresses  through  the  convolute  and  the  uniform 
transmission  line.  Each  feature  will  be  discussed 
and  then  the  results  of  the  Hite  calculations  will 
be  presented. 


where  Ym  is  the  position  of  the  sheath,  YQ  is  the 
position  of  the  anode  and  the  cathode  is  at  Y  ’  0. 
The  values  of  V  at  four  positions  equally  spaced 
about  (x,y)  are  calculated  and  the  electric  field 
is  calculated  from 

E  -  -VV  (10) 

The  local  magnetic  field  B  =  Bz  z  is  given  by 

me2 

Bz  =  -—  (Y-  1)  0<Y<Ym  (11a) 


The  test  particles  are  assumed  to  originate  in  a 
cathode  plasma  at  a  y  coordinate  yQ  -  10~&  m, 
where  the  voltage  is  <_  0.3  eV  and  the  magnetic 
vector  potential  is  AQ ,  which  is  calculated  from 
the  parapotential  theory.2  The  initial  energy  of 
the  electron  is  chosen  between  1  and  10  eV  to 
simulate  the  effect  of  electron  emission  from  a 
few  eV  plasma.  The  initial  particle  energy  deter¬ 
mines  che  absolute  value  of  the  electrons  initial 
v_elocit^_  UQ.  The  initial  canonical  momentum 
P  =  P.c  x  is  assumed  to  be  zero,  so  initially. 


u 

x  7  m 
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(3) 


Bz  - 


11  oh 
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if  Y  <  Y  <  Y 
m  —o 


(lib) 


The  relativistic  equation  of  motion  for  the 
electrons  is  j  _  _ 

(73HU)  «  -e  (E  +  IT  X  B) 


(12) 


and  is  combined  with  the  local  electric  and 
magnetic  fields.  It  is  then  solved  using  the 
integration  routine  STF0DE°  to  find  the  new 
particle  postion  and  velocity  components,  within  a 
given  error  criteria,  after  a  time  increment. 

The  particle  is  progressively  accelerated  through 
the  transition  and  transport  sections  of  line  for 
successive  time  steps. 


Results 

The  Mite  experiment2’^  used  two  transition 
geometries  to  change  the  gap  spacing  from  0.02  to 
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0.01  m  in  a  triplate  vacuum  transmission  line  with 
an  effective  width  w  «  0.50  m.  In  the  first 
geometry,  the  transition  was  made  over  an  axial 
length  of  0.04  m  and  40  percent  of  the  power  was 
lost  between  0.50  and  1.4  m  from  the  beginning  of 
the  uniform  line.  The  second  transition  was  made 
over  an  axial  length  of  0.14  m  and  the  power  trans¬ 
port  was  about  100  percent  efficient.  These  two 
geometries  were  simulated  with  typical  Mite  para¬ 
meters  of  VQ  =  2.0  MV,  IT  =  0.4  MA. 

The  effect  of  having  space  charge  in  the  vacuum  gap 
is  illustrated  in  Fig.  2.  The  equipotentials  for 
a  1  cm  taper  are  shown  for  the  Mite  parameters  of 
V  and  1^.  The  position  of  the  edge  of  the 
electron  sheath  is  Y_  and  is  shown.  The  effect 

IQ 

of  the  space  charge  is  to  produce  a  positive  Ex 
near  the  cathode  and  to  distort  the  distribution 
of  Ey  and  B, . 


Fig.  2.  Effect  of  space  charge  in  vacuum  gap  of 

1  cm  long  transition  is  shown.  The  dotted 
lines  are  equipotential  at  200  kV  intervals 
and  Ym  is  the  electron  sheath  position. 

For  the  severe  0.01  m  long  transition  the  distor¬ 
tion  is  not  very  large.  The  maximum  value  of 
EX/(V  /d)  is  only  0.005  for  the  1  cm  taper  and  is 
much  less  for  the  0.04  m  and  0.14  m  tapers.  Con¬ 
sequently,  the  effect  of  the  convolute  on  the 
electron  motion  is  small  and  must  be  calculated 
with  a  very  small  relative  and  absolute  errors  of 
8  =  10“6  X. 


The  emission  current  as  a  function  of  X  is  also 

0 

shown.  Since  the  differential  electron  current 
Jg(x)  =AIg/AX  is  supplied  from  the  cathode,  the 
approximate  shape  of  the  canonical  momentum  dis¬ 
tribution  F(PX)  can  be  approximated  by 

1  ^?x ^ 

F(PX)~  y  wee  dP  /dx 
m  x 

which  is  shown  in  Fig.  6  for  both  tapers  under  the 
assumption  that  the  initial  energy  W  of  the  elec¬ 
trons  is  10  eV  at  the  cathode  surface. 


Fig.  3.  The  electron  current  I  and  the  calculated 
emission  current  per  unit  length  for  the 
0.14  m  long  transition  is  shown. 
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Fig.  4.  The  final  canonical  momentum  Px  and  the 
electron  emission  current  Jg  vs.  the 
initial  electron  position  XQ  for  the 
0.04  m  taper. 


The  current  Ig  carried  by  the  electron  flow  increase 
steadily  as  the  spacing  between  conductors  is 
reduced  in  the  transition  section,  as  shown  in 
Fig.  3  for  the  0.14  m  taper.  The  final  canonical 
momentum  that  the  electron  achieves  as  it  is 
accelerated  through  the  convolute  is  shown  as  a 
function  of  its  initial  position  XQ  in  Fig.  4  and 
5  for  the  0.04  m  and  0.14  m  convolutes  respectively. 


Discussion  and  Conclusions 

In  both  cases,  the  canonical  momentum  is  negative 
and  F(PX)  has  a  very  small  width.  The  spread  in 
Px  is  %10”°  of  that  allowed  in  the  uniform  line. 
The  small  values  of  Px  obtained  with  space  charge 
are  much  less  than  the  values  estimated  from  the 
vacuum  fields  alone. 
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Fig.  8.  The  calculated  distributions  F(P  )  for  the 
Mite  transition  sections. 


The  fact  that  ?  <  0  for  the  injected  electrons 

means  that  they  will  flow  for  many  Larmor  radii. 
Electron  with  W  =  10  eV  travel  through  the  line 
past  the  0.50  m  position  at  which  the  losses  occur 
and  could  be  susceptible  to  an  instability  with  a 
spatial  growth  length  much  less  than  0.50  m,  as 
hypothesized4  to  explain  the  Mite  results. 

Since  the  electron  current  is  a  very  strong  func¬ 
tion  of  the  separation  between  plates,  the  bulk  of 
the  electrons  in  the  flow  have  small  values  of  Px. 
Even  though  the  gradients  in  the  fields  are  larger 
near  the  end  of  the  transition,  the  electrons  are 
accelerated  through  the  transition  section  before 
they  acquire  a  large  negative  canonical  momentum. 


The  electrons  with  more  negative  values  of  Px 
originate  further  from  the  output  of  the  transition 
section.  The  effect  of  using  a  more  gradually 
tapered  transition  section  is  to  broaden  the  canon¬ 
ical  momentum  distribution.  Since  the  more  gradual 
taper  has  efficient  power  propagation,  the  results 
indicate  that  a  broader  F(PX)  provide  more  reliable 
transport  in  long  lines.  An  injector  designed 
such  that  Jx  approximately  equals  a  constant 
through  a  long  transition  section  should  be  the 
optimum  arrangement.  Such  a  transition  section 
will  be  designed  and  tested  on  the  Mite  experiment. 

In  conclusion,  the  simulations  indicate  that  the 
difference  between  the  distributions  F(PX)  for  the 
lossy  and  the  efficient  transitions  is  small  but 
significant.  The  results  suggest  a  way  to  improve 
the  transition  and  define  experiment  and  theory 
development  required  to  explore  the  implications 
further.  The  results  indicate  that  an  experiment7 
to  measure  F(PX)  for  the  Mite  transition  sections 
should  be  capable  of  resolving  the  distributions 
in  Fig.  6.  Finally,  the  stability  of  electron 
flow  with  F(P  )  similiar  to  those  presented  in 
Fig.  6  should  be  examined  under  the  conditions  of 
the  large  E  and  Bz  present  in  self  magnetically 
insulated  transmission  lines . 
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